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Abstract 

The production of transverse energy in deep inelastic scattering is measured as a 
function of the kinematic variables x and Q 2 using the HI detector at the ep collider 
HERA. The results are compared to the different predictions based upon two alter- 
native QCD evolution equations, namely the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) and the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equations. In 
a pseudorapidity interval which is central in the hadronic centre of mass system 
between the current and the proton remnant fragmentation region the produced 
transverse energy increases with decreasing x for constant Q 2 . Such a behaviour 
can be explained with a QCD calculation based upon the BFKL ansatz. The rate 
of forward jets, proposed clS cl SI gnature for BFKL dynamics, has been measured. 
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1 Introduction 



The electron-proton collider HERA has opened new kinematical regions in the study 
of Deep Inelastic Scattering (DIS): the regions of large four-momentum transfer Q 2 (up 
to Q 2 ~ 10 4 GeV 2 ) and small Bjorken-x (down to x ~ 10~ 4 ). It has been suggested 
that the small x region may be sensitive to new dynamic features of QCD, e.g. 0. The 
ZEUS and HI collaborations have observed 0, [J that the proton structure function 
F 2 exhibits a strong rise towards small Bjorken-x. This rise has caused much debate 
on whether the HERA data are still in a regime where the QCD evolution of the parton 
densities can be described by the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) 0] 
evolution equations, or whether they extend into a new regime where the QCD dynamics 
is described by the BFKL (Balitsky-Fadin-Kuraev-Lipatov) || evolution equation. The 
BFKL evolution equation is expected to become applicable in the small x region, since 
it resums all leading a s In 1/x terms in the perturbative expansion, in contrast to the 
DGLAP equation. Present F 2 measurements do not yet allow to discriminate between 
BFKL and conventional DGLAP dynamics || [7|], and are perhaps too inclusive a measure 
to be a sensitive discriminator. Hadronic final states may give additional information and 
could be more sensitive to the parton evolution || || pL0| . In this paper we use the HI 
detector to study properties of the hadronic final state, namely the transverse energy 
flow and the production of forward jets, to test the BFKL versus DGLAP hypothesis. 
Throughout this paper, "forward" refers to the region around the proton direction in 
the HERA laboratory frame, and transverse quantities are calculated with respect to the 
proton direction. 

For events at low x, hadron production in the region between the current jet and 
the proton remnant is expected to be sensitive to the effects of the BFKL or DGLAP 
dynamics. At lowest order the BFKL and DGLAP evolution equations effectively resum 
the leading logarithmic a s In 1/x or a s In Q 2 contributions respectively. In an axial gauge 
this amounts to a resummation of ladder diagrams of the type shown in Fig. |l]. This 
shows that before a quark is struck by the virtual photon, a cascade of partons may be 
emitted. The fraction of the proton momentum carried by the emitted partons, Xi, and 
their transverse momenta, k^i, are indicated in the figure. In the leading log DGLAP 
scheme this parton cascade follows a strong ordering in transverse momentum k\ n ^> 
^Tn-i 3> ••■ 3> ^ti; while there is only a soft (kinematical) ordering for the fractional 
momentum x n < x n -i < ... < x\. In the BFKL scheme the cascade follows a strong 
ordering in fractional momentum x n -C £ n -i ^ ••• Xi, while there is no ordering in 



transverse momentum JO]]. The transverse momentum follows a kind of random walk in 
k,T space: the value of k,Ti is close to that of /cti-i, but it can be both larger or smaller Jl2| . 
As a consequence, BFKL evolution is expected to produce more transverse energy Et than 
DGLAP evolution [JT^, |14]] in the region between the struck quark and the remnant for low 
x events. Intriguingly, a numerical BFKL calculation of the transverse energy flow [|L~f 
came out close to the measurements published by HI [T5|] , while the DGLAP prediction 
is too low. In addition, the DGLAP calculation predicts an increase in transverse energy 



with rising x, while the BFKL calculation predicts the opposite |13|. With increased 



statistics, the x dependence of this effect can be studied in more detail in this paper. 
Another possible signature of the BFKL dynamics is the rate of jets with transverse 
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Ej et /E p , large com- 



Here E ]et and E p are the energies of the jet and the 



momentum kxjet ~ Q and the momentum fraction of the jet, Xj et 
pared with Bjorken-x || [1C], |16| 
incoming proton respectively. Due to the strong ordering in DGLAP evolution, the con- 
dition A/rjet 2 ps Q 2 suppresses the phase space for jet production. However jet production 
from BFKL evolution is governed by the ratio Xj et /x, which is large. Hence the rate of 
events with a jet satisfying the selection is predicted to be higher for the BFKL than for 
the DGLAP scenario. 

Apart from numerical calculations, predictions for final state observables are also avail- 
able as Monte Carlo models, based upon QCD phenomenology. In this report we consider 
two of the currently available Monte Carlo programs: the MEPS (Matrix Elements plus 
Parton Showers) and CDM (Colour Dipole Model) models. Both provide good descrip- 
tions of fixed target DIS and e + e~data [[17]]. The CDM model |18| provides an implemen- 
tation of the colour dipole model of a chain of independently radiating dipoles formed by 
emitted gluons |19| . Photon-gluon fusion events are not described by this picture and are 
added at a rate given by the QCD matrix elements [2C|. The CDM description of gluon 
emission is similar to that of the BFKL evolution, because the gluons emitted by the 
dipoles do not obey strong ordering in kx ||21|| . The CDM does not explicitly make use 
of the BFKL evolution equation, however. The MEPS model is an option of the LEPTO 



generator |2(J based on DGLAP dynamics. MEPS incorporates the QCD matrix elements 
up to first order, with additional soft emissions generated by adding leading log parton 
showers. The emitted partons are strongly ordered in kx- Both Monte Carlo programs use 



the Lund string model [22] for hadronizing the partonic final state. The parton density 
parametrization used here is that of MRSH | 23fl , which results in a good description of 
the HERA F 2 measurements H, HI. 





X P ^Tl~Q 



jet 



Figure 1: Parton evolution in the ladder approximation. The selection of forward jets in DIS events is 
illustrated. 
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2 Detector Description 



A detailed description of the HI apparatus can be found elsewhere ||24j| . The following 
briefly describes the components of the detector relevant to this analysis. 

The hadronic energy flow and the scattered electrons are measured with a liquid 
argon (LAr) calorimeter and a backward electromagnetic lead-scintillator calorimeter 
(BEMC). The LAr calorimeter [2~5] extends over the polar angular range 4° < 9 < 153° 
with full azimuthal coverage, where 9 is defined with respect to the proton beam direction 
(+z axis). It consists of an electromagnetic section with lead absorbers and a hadronic 
section with steel absorbers. Both sections are highly segmented in the transverse and 
longitudinal directions with about 44000 cells in total. The total depth of both calorime- 
ters varies between 4.5 and 8 interaction lengths for 9 < 125°. Test beam measurements of 



the LAr calorimeter modules show an energy resolution oIoeIE « 0.50/ yE [GeV] ©0.02 
for charged pions [p6fl . The hadronic energy measurement is performed by applying a 
weighting technique in order to account for the non-compensating nature of the calorime- 
ter. The absolute scale of the hadronic energy measurement is presently known to 5%, as 
determined from studies of the transverse momentum (pt) balance in DIS events. 

The BEMC (depth of 22.5 radiation lengths or 1 interaction length) covers the back- 
ward region of the detector, 151° < 9 < 176°. The BEMC energy scale for electrons 
is known to an accuracy of 1.7%. Its resolution is given by <Je/E = 0.10/ yE [GeV] © 



0.42/£[GeV] ©0.03 |27 



The calorimeters are surrounded by a superconducting solenoid which provides a uni- 
form magnetic field of 1.15 T parallel to the beam axis in the tracking region. Charged 
particle tracks are measured in a central drift chamber and the forward tracking system, 
covering the polar angular range 7° < 9 < 165°. A backward proportional chamber 
(BPC), in front of the BEMC with an angular acceptance of 155.5° < 9 < 174.5° serves 
to identify electrons and to precisely measure their direction. Using information from 
the BPC, the BEMC and the reconstructed event vertex the polar angle of the scattered 
electron is known to a precision of 2 mrad. 



3 Event Selection and Kinematics 

The data used in this analysis were collected in 1993, with electrons of energy E e = 
26.7 GeV colliding with protons of energy E p = 820 GeV, resulting in a total centre of 
mass energy of y^s = 296 GeV . The data correspond to an integrated luminosity of 320 
nb _1 . For this analysis DIS events with Q 2 < 100 GeV 2 are used, in which the scattered 
electron is observed in the BEMC. The events are triggered by requiring a cluster of more 
than 4 GeV in the BEMC. After reconstruction, DIS events are selected in the following 
way: 

• The scattered electron, defined as the most energetic BEMC cluster, must have an 
energy E' e larger than 12 GeV and a polar angle 9 e below 173° in order to ensure 
high trigger efficiency and a small photoproduction background || . 
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The lateral size of the electron cluster, calculated as the energy weighted radial 
distance of the cells from the cluster centre, has to be smaller than 4 cm. The 
cluster must be associated with a reconstructed BPC space point which must lie 
within 4 cm of the cluster centre of gravity. Further reduction of photoproduction 
background and the removal of events in which an energetic photon is radiated off 
the incoming electron is achieved by requiring J2j (Ej — p z j) > 30 GeV [[|, with 
the sum extending over all particles j (measured calorimetrically) of the event. 

The radial coordinate of the BPC hit must be less than 60 cm, corresponding to an 
electron angle above 157° with respect to the nominal interaction point, ensuring 
full containment of the electron shower in the BEMC. 

The z position of the event vertex reconstructed from charged tracks has to be 
within 30 cm of the average of all collision events. 

The energy in the forward region (4.4° < 6 < 15°) has to be larger than 0.5 GeV 
in order to exclude diffract ive- like events with large rapidity gaps in the forward 
region 



28, 15 . 



The kinematic variables are determined using information from the scattered electron: 
Q 2 = 4E e E' e cos 2 (6 e /2) and y = 1 - (E'jE e ) ■ sin 2 (fl e /2). The scaling variable x is then 
derived via x = Q 2 / (ys), and the hadronic invariant mass squared is W 2 = m 2 + sy — Q 



2 



As the precision of the y measurement degrades with 1/y, a cut y > 0.05 is imposed. 
Events suffering from QED radiation or from a badly reconstructed electron are 
removed by requiring that they also fulfil this cut if y is calculated from the measured 
hadrons. 



4 Transverse Energy Flows 

The event sample for the energy flow measurements, in which 60% of the total luminosity 
has been used, contains 9529 DIS events with 5 GeV 2 < Q 2 < 100 GeV 2 and with 
10~ 4 < x < 10~ 2 . The measurements are performed in the hadronic centre of mass 
system (CMS), that is the rest system of the proton and the exchanged boson. The 
orientation of the CMS is such that the direction of the exchanged boson defines the 
positive z' axis. 

The x - Q 2 plane is divided into three slices with constant Q 2 and varying x, re- 
sulting in 9 kinematic bins (see Table [j]). In Fig. ^|, the flow of transverse energy, as 
measured with the LAr and BEMC calorimeters, is shown as a function of pseudorapidity 
i] = — lntan(#/2) for the individual kinematic bins. The data are corrected for detector 
effects using the CDM generator and a full simulation of the HI detector response. The 
data in general exhibit a moderate slope, extending from the current region towards the 
proton remnant up to the edge of the detector acceptance. For large x and Q 2 the data 
are reasonably well described by the two models under consideration, MEPS and CDM. 
At small x and Q 2 however, the models, in particular MEPS, deviate significantly from 
the data. The MEPS model generates too little Et- The CDM is able to describe the 
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2.12 ±0.12 
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1.99 ±0.13 




3.0-10 


20-50 


4.9 


32.6 


1.75 ±0.17 





Table 1: Values of (Et ) measured in the CMS pseudorapidity interval from —0.5 to 0.5 as a function 
of x and Q 2 . The errors quoted contain the statistical and the systematic point-to-point errors added in 
quadrature. An overall uncertainty of 9% is not included. Also given are perturbative BFKL predictions 
calculated according to |Q (see text). 

Et flow towards the remnant, but overestimates the Et in the current region, around 
7] = 3. We note that in events where the hard subprocess leads to two hard jets in the 
detector, both models provide a good description of the energy flow [^9j around the jets. 

In order to quantify the evolution of Et with x and Q 2 , the average Et per unit of 
pseudorapidity is measured in the central pseudorapidity range —0.5 < r\ < 0.5. Note that 
this corresponds roughly to the forward range 2 < rji ah < 3 in the laboratory frame. The 
measured values of (Et ) as a function of x are shown for the three Q 2 slices in Fig. ^a, 
and summarized in Table |I[ They are of the order of 2 GeV per unit of pseudorapidity. 
For (Q 2 ) w 14 GeV 2 , the (E T ) drops by about 25% when going from (x ) = 0.37- 10~ 3 to 
2.3-10~ 3 . Though less significant, the data at (Q 2 ) w 8 GeV 2 and at (Q 2 ) w 30 GeV 2 
confirm these (Et ) measurements. For fixed x, the level of Et rises slightly with Q 2 Q. 

The errors shown in Fig. [|a are the statistical and systematic point-to-point errors 
added in quadrature, which are typically 5 % each. The corrections for detector effects are 
typically 20%. The uncertainty in the BEMC calibration (1.7%) leads to a sizeable effect 
(4%) in the bin of largest x and Q 2 . Small effects from QED radiation (typically 1-2%) 
and from photoproduction background in the lowest x and Q 2 bin (3%), estimated from 
Monte Carlo simulations, are not corrected for but absorbed in the systematic error. The 
model dependence of the correction was investigated with the CDM and MEPS models 
and leads to a 4% point-to-point error and a 5% overall error. In addition, an overall 
scale error of 7% arises from the calorimeter calibration (5%), and from varying details 
of the analysis method, such as the treatment of noise in the calorimeter, clustering of 
calorimeter cells, and the simulation of the calorimeter response, affecting the result by 
5% in the forward region. 

Fig. |3|a and Table [I] display the prediction for (Et ) from a BFKL based QCD calcu- 

1 When W and Q 2 are chosen as kinematic variables, (Et ) is found to rise with W and to be almost 
independent of Q 2 . 
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Figure 2: The flow of transverse energy Et in the CMS as a function of pseudorapidity -q. The 
remnant direction is to the left. Shown are data for nine different kinematic bins (see Table 1). The 
distributions are normalized to the number of events N in each case. Only statistical errors are shown. 
For comparison, the models CDM (full line) and MEPS (dashed) are overlayed. 
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lation |TJ|], yielding values of the order of 1.7 GeV per unit of pseudorapidity and with 



similar slopes as a function of x as seen in the data. The validity of the calculation [J4 
is limited to the range shown. A DGLAP based calculation yields an (Et ) of around 
0.4 GeV jnj with the opposite x dependence. These calculations however do not include 
the non-perturbative hadronization. This can presently only be included using Monte 
Carlo models which contain both the perturbative QCD evolution and a phenomenolog- 
ical hadronization model. Since the CDM provides a reasonable description of the data, 
and since the perturbative BFKL calculation of (Et ) is close to the partonic final state 
(Et ) produced perturbatively in the CDM (see Fig. |]a), the CDM may provide an es- 
timate of the effect of hadronization on the BFKL (Et ) calculation. That effect can 
be seen in Fig. ^]a, where the (Et ) according to the CDM is plotted before and after 
hadronization. It amounts to an increase in (Et ) of 0.3-0.4 GeV and is rather indepen- 
dent of x. When taking into account the hadronization effect as modelled by the CDM, 
the BFKL calculation is in good agreement with the data. In this context it is interest- 
ing to observe that the Monte Carlo model without fc^ ordering, the CDM, is consistent 
with the BFKL calculation for the perturbative part, and that it is able to describe the 
measured magnitude and x dependence of (Et ) reasonably well. 

It has to be stressed that the model predictions depend on a variety of parameters 
which can be tuned fl5| . For this analysis the model settings were chosen such as to render 



a good description of the measured energy flow at large x and Q 2 , a kinematic region 
where theoretical uncertainties are minimal. CDM and MEPS were used in the program 
implementations ARIADNE 4.03 and LEPTO 6.1, respectively. In the MEPS model, 
divergences of the matrix element are avoided with a cut-off for parton-parton invariant 
masses, m^- > y cnt ■ W . For this analysis the MEPS cut-off was parametrized such as to 
follow the limit at which the order a s contribution exceeds the total cross section within 
a margin of 2 GeV . For the CDM the standard value of the parameter y cut (y cu t = 
0.015) is used to regulate the admixture of boson-gluon fusion events according to the 
matrix element. Turning off the boson-gluon fusion admixture, or using the same cut-off 
prescription as described above has little influence on the CDM predictions. 

Improvements to the MEPS model are conceivable, for example with an improved 
scheme for matching parton showers and matrix elements, by changing the arrangement 
of color connections, or by changes to the remnant fragmentation which is not well tested 
at small x. Whether or not such improvements can result in an acceptable description 
of the energy flows while maintaining a good description of other final state observables 
can presently not be judged. Therefore the failure of the MEPS model to describe (Et ) 
cannot be unambiguously identified with the fact that its parton shower evolution is based 
upon the DGLAP evolution. For the same reasons the success of the CDM in this respect, 
where the parton cascade does not obey strong /ct ordering as in the BFKL evolution, is 
intriguing but may be fortuitous. 

Having identified the average Et as a sensitive observable to test QCD predictions, it 
is interesting as well to study the shape of the Et distribution. The distribution of the 
observed Et in the central pseudorapidity bin —0.5 < i] < 0.5 is shown in Fig. ^]b for a 
kinematic bin at low x and Q 2 . The CDM gives a reasonable description of the data, and 
the MEPS model produces events accumulating at lower values of Et than the majority 
of the data events. The difference in (Et ) between the data and the MEPS model does 
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Figure 3: a) Transverse energy as a function of x for three different values of Q 2 . The transverse energy 
(Et ) is measured in the CMS in the pseudorapidity range from —0.5 to +0.5. Shown are the HI data, 
corrected for detector effects, the BFKL calculation and the prediction from the CDM before (partons) 
and after (hadrons) hadronization. The error bars contain the statistical and the systematic point-to- 
point errors added in quadrature. An overall scale uncertainty of 9% is not shown, b) The uncorrected 
distribution of the Et (as defined above) for the events of the kinematic bin with (x ) = 0.37T0 -3 
and (Q 2 ) = 13.1 GeV 2 . The data are compared with the CDM and MEPS models, including a full 
simulation of the HI detector. 



not seem to stem from a tail in the data. 



5 Forward Jets 



We have studied DIS events at small x which have a jet with large Xj C t ||30|| . A cone 
algorithm is used to find jets, requiring an Et larger than 5 GeV in a cone of radius 
R = y/Arj' 2 + A0 2 = 1.0 in the space of pseudo-rapidity r] and azimuthal angle <fi in the 
HERA frame of reference. In addition to the selections given in section 3, the requirement 
y > 0.1 was imposed to ensure that the jet of the struck quark is well within the central 
region of the detector and is expected to have a jet angle larger than 60°. In this sample 
of DIS events with Q 2 w 20 GeV 2 and 2-10 _4 <x<2-10~ 3 we have counted events which 
have a "forward" jet defined by x iet > 0.025, 0.5 < p T ]et 2 /Q 2 < 4, 6° < jot < 20° and 
PTjet > 5 GeV, where prjet is the transverse momentum of the jet. A typical event with 
a high energy forward jet is shown in Fig. f|a. The transverse energy flow around the 
forward jet axis, averaged over all selected events, is shown versus r\ and 4> in Figs. Hb 
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x range 


data 
events 


MEPS 

events 


CDM 

events 


<j(ep — > jet + X) 
(pb) 


2 • 10" 4 - 1 • 1(T 3 


271 


141 


282 


709 ± 42 ± 166 


1 • io~ 3 - 2 ■ 10~ 3 


158 


101 


108 


475 ±39 ±110 



Table 2: Numbers of observed DIS events with a selected forward jet, corrected for radiative events 
faking this signature. These may be directly compared with the expectations from the Monte Carlo 
models. The measured cross section ep — > jet + X for forward jets is also given. The errors reflect the 
statistical and systematic uncertainties. 



and |]c. Distinct jet profiles are observed, which are well described by the Monte Carlo 
models. 

The resulting number of events observed with at least one forward jet in the kine- 
matical region 160° < 9 e < 173° and E e > 12 GeV is given in Table Q and compared 
to expectations of the MEPS and CDM models after detector simulation. The data are 
corrected for photoproduction background and radiative events, which due to the changed 
kinematics at the hadron vertex can eject a jet in the forward direction. About 4% of 
the data events were found to contain two forward jets. Results on jet production for 
both models have been compared with data in previous analyses [[!l], |32|]. In particular 



the MEPS model was found to show good agreement with the data outside the forward 
region. In the kinematic range studied here the CDM generally describes the data better 
than the MEPS model. The predictions were found not to depend significantly on the 
uncertainties of the proton structure function. However, increasing the Xj e t cut from 0.025 
to 0.05 reduces the total number of events with forward jets to 46 for CDM, to 77 for 
MEPS and to 105 for the data, hence CDM does not describe the rate of high energy jets. 

The measured cross section for forward jets satisfying the cuts given above is also 
presented in Table |2|. It has been corrected for detector effects using the CDM. The sys- 
tematic errors include effects from DIS event selection, the calorimeter energy scale (5%), 
the jet angle bias (10 mrad), the proton structure function, and a global normalization 
uncertainty of 4.5%. Event pile-up effects were found to be negligible. The systematic 
errors on the two data points are largely correlated. The ratio of the jet cross section for 
the low x to the high x bin is 1.49 ± 0.25. 

The precision of the data does not yet allow a firm conclusion to be drawn. We note, 
however, that the forward jet cross section is larger in the low x bin than in the high x bin. 
This is expected from BFKL dynamics as an analytical calculation |T0, |33| at the parton 
level demonstrates: in the kinematical region selected the ratio of the cross sections in 
the low x bin to high x bin is 1.62 for a calculation including BFKL evolution, compared 
to 1.03 for a calculation without gluon emission from the ladder in Fig. [l|. 



6 Conclusions 

In order to shed light on the QCD mechanism responsible for parton evolution in the 
regime of small Bjorken-x, the production of particles and of jets in the forward region 
has been measured at HERA. A forward jet selection designed to enhance the yield in 
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Figure 4: a) DIS event event with a forward jet in the HI detector. The protons are incident from 
the right, electrons from the left. The scattered electron is detected in the backward electromagnetic 
calorimeter (BEMC) with an angle of 166° and an energy of 18.9 GeV. The forward jet is observed in 
the liquid argon calorimeter and has an angle 6*j e t = 11° and an energy £j e t = 65 GeV. Averaged over all 
events with a selected forward jet, the transverse energy flow around the forward jet axis is shown in b) 
as a function of Ar), integrated over \ A(f>\ < 1.0 and in c) as a function of A</>, integrated over | Ary| < 1.0. 
Here A-q and A</> are measured with respect to the reconstructed jet axis. 
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the case of BFKL evolution, and to suppress the yield for DGLAP evolution, results in a 
rate of observed forward jets compatible with the BFKL expectation. A firm conclusion 
on the growth with x however would necessitate a larger data sample. 

The flow of transverse energy versus pseudorapity in the hadronic centre of mass 
system has been measured as a function of the kinematic variables x and Q 2 . The ob- 
served magnitude and x dependence of the average transverse energy is in agreement 
with a perturbative QCD calculation based upon the BFKL mechanism, assuming a non- 
perturbative contribution from hadronization as predicted by the colour dipole model. 
Though the BFKL mechanism provides a natural explanation of the data, it is also pos- 
sible that the currently used hadronization scheme is inadequate, and that improved 
DGLAP based models may also be able to describe the observed features of transverse 
energy production. The data, which are corrected for detector effects, provide important 
constraints for further development of the understanding of QCD at small x. 
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